Non-Newtonian gravity in finite nuclei 
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Abstract. In this talk, we report our recent study of constraining the non-Newtonian gravity 
at femtometer scale. We incorporate the Yukawa-type non-Newtonian gravitational potential 
consistently to the Skyrme functional form using the exact treatment for the direct contribution 
and density-matrix expansion method for the exchange contribution. The effects from the non- 
Newtonian potential on finite nuclei properties are then studied together with a well-tested 
Skyrme force. Assuming that the framework without non-Newtonian gravity can explain the 
binding energies and charge radii of medium to heavy nuclei within 2% error, we set an upper 
limit for the strength of the non-Newtonian gravitational potential at femtometer scale. 



1. Introduction 

How to understand the nature of the gravitational force and unify it with other fundamental 
forces are among the most important questions to answer in the new century pQ. It has gradually 
been realized that the traditional Newtonian gravitational potential, which obeys the Inverse- 
Square-Law (ISL), can not explain all the phenomena relative to gravity in nature. Great efforts 
have been devoted to search for the possible existence of non-Newtonian gravity that violates 



the ISL [21 [31 HI [5j El El [8] . Especially, these studies have put upper limits of the strength of 
the non-Newtonian gravity from as large as galaxy scale to as small as 10 fm [7]. At galaxy 
scale, various modified gravity theories are used to study the non-Newtonian gravity such as 
the scalar-tensor-vector gravity and f(R) gravity. At short distances, the origin of the non- 
Newtonian potential can be attributed to the possible existence of extra dimensions within 
string/M theories or exchanging light and weakly-coupled new particles in the supersymmetric 
extension of the Standard Model. 

The following form of the gravitational potential between two objects of masses mi and mi 
at positions r\ and r-i has been widely used at both galaxy scale and short distances [9] 

Wrl,r 2 ) = + (1) 

In - r 2 | 

In the above, the first term is the traditional Newtonian potential where G is gravitational 
constant, while the second term represents the non-Newtonian potential, where a and A are the 
strength and length scale parameters, respectively. 

At galaxy scale (A ~ 10 21 m), this simple Yukawa- type non- Newtonian gravitational potential 
can be viewed as a reduced form of the scalar-tensor- vector gravity or the f(R) gravity in the 
weak-field limit [101 [TT] , and it has been used to successfully explain both the flattening of the 
galaxy rotation curves away from the Kepler limit |12l [T3] and the Bullet Cluster 1E0657-558 
observations in the absence of dark matter |14j . For a recent review of applying this Yukawa- 
type non-Newtonian gravitational potential on well-known observations at galaxy scale, we refer 
the readers to Ref. [T5] , 

At short distances, the Yukawa potential in Eq. ([1]) may come from exchanging a light 
and weakly-coupled spin-0 axion [16] or spin-1 U-boson [T7] corresponding respectively to an 
attractive or repulsive potential. The exchanged boson is related to several interesting new 
phenomena in particle physics and cosmology |18|. [19] and may mediate the annihilation of dark 
matter particles [20, 21] which accounts for the 511 keV 7-ray from the galactic bulge as well. 

Especially, the properties of a neutron star would be affected by the possible existence of 
the non-Newtonian gravitational potential if its length scale is smaller than the size of the 
star. In this case, if one still wants to rely on solving the Tolman-Oppenheimer-Volkoff using 
the equation of state (EOS) of neutron star matter, the above Yukawa-type non-Newtonian 
gravitational potential can be viewed as an extra potential between nucleons in addition to 
the basic nucleon-nucleon interaction. It has been found that the possible existence of the 
repulsive Yukawa potential from exchanging a spin-1 U-boson would stiffen the EOS of neutron 
star matter compared to that from a nucleon-nucleon interaction only and thus increase both 
the mass and the radius of a neutron star |22j . Even if with a supersoft symmetry energy at 
suprasaturation densities, the repulsive Yukawa potential could still maintain a stable neutron 
star [23]. In addition, the core-crust transition density /pressure could also be affected and the 
neutron star structure could be modified in the presence of the Yukawa potential [23]. The 
coupling of U-boson and nucleon has been successfully introduced into the relativistic mean- 
field model and similar effects have been observed |25] , Furthermore, it is shown that although 
the presence of hyperons and quarks would soften the EOS of neutron star matter, the repulsive 
Yukawa potential can stiffen the EOS and account for the heavy mass of PSR J1614-2230 |26j. 
Thus, our knowledge of the non-Newtonian potential and the nuclear EOS are both important 
in understanding many interesting phenomena in nuclear astrophysics. 

At a certain length scale A, the non-Newtonian potential is generally compared with the 
typical interaction for the corresponding length scale, through which the strength of the non- 
Newtonian potential is constrained if the interaction is well determined. As mentioned above, 
the study and the constraint of the non-Newtonian potential below 10 fm are still lacking. The 
most typical system at femtometer scale is a finite nucleus, where the typical interaction is the 



strong interaction. Thanks to the great efforts made by the nuclear physicists, the properties 
of finite nuclei can be very well explained by effective nucleon-nucleon interactions. Especially, 
from the empirical values of isoscalar and isovector macroscopic quantities, the parameters in the 
Skyrme effective interaction can be inversely determined, leading to the MSLO Skyrme force |27j . 
This gives us the opportunity to study the possible existence of the Yukawa-type non-Newtonian 
potential in the finite nuclei whose properties have been well described by the MSLO Skyrme 
force, and constrain the strength of the non-Newtonian potential at femtometer scale. 



2. Incorporating the Yukawa potential to the Skyrme-Hartree-Fock calculation 

In the one-boson-exchange picture, the Yukawa-type non-Newtonian gravitational potential 
between two nucleons at position r\ and r*2 is written as 

V Y (n,f 2 ) = ±^ (2) 

47r |n — r2\ 

In the above, g is the boson-nucleon coupling constant and it can been expressed as g = 
\J 4Tr\a\Gm 2 comparing Eq. ([2]) with Eq. (pQ), where m is the nucleon mass, p, = 1/A is the 
mass of axion or U-boson which is largely quite uncertain, and it determines the length scale of 
the non-Newtonian potential. The potential in Eq. ([2]) can be positive or negative, depending on 
the spin of the exchanged boson. To ease discuss we omit the ± sign in the following formulas. 

To apply the above Yukawa potential to finite nuclei calculation, we need to do Hartree-Fock 
calculation as for the Skyrme interaction. The additional potential energy due to the existence 
of the Yukawa potential can thus be written as 

E Y = \ X>J IW - PrP*P T )\ij), (3) 

where P r , P a , and P T are the space, spin, and isospin exchange operator, respectively, and \i) is 
the quantum state of the ith. particle containing spatial, spin, and isospin parts. The coefficient 
^ is due to the double counting for exchanging i and j. 

The first term in Eq. ([3]) is the direct (Hartree) contribution. By expressing the quantum 
state in the spatial coordinate representation it is written as 

1 r a 2 p~ lA^i— ^2 1 



Ey = » / Pin)p(r 2 )^-— ^cfVicfra, (4) 

2 J 4-7T \r\ — V2\ 

where 

p(r) = Y,<t>Ti(?,°)<f>Ti(r,cr) (5) 

i,(T,T 

is the nucleon number density with <f) T i(f,a) being the spatial wave function of the ith. particle 
with spin a and isospin r. 

The second term in Eq. is the exchange (Fock) contribution. By using P a = (l + c?i -(?2)/2 
where <?i(2) is the Pauli operator acting on the state \i) it can be written as 

E Y = -7 Y] / [Pr(rL,r 2 )/9 T (r 2 ,fi) + p" T (fi,f 2 ) ■ p" T (r2,ri)]- p ^—d 3 rid 3 r2, (6) 



r=n,p 



where 



/5 r (n,f 2 ) =Y4>Ti{n,<?)4>Ti{r2,<?) (7) 



1,(7 



and 

Pr{ri,r 2 ) = J2J2(j) i ; i {f 1 ,a')(t) T i{f2,a-)(a'\a\a) (8) 

i a, a' 

are the off-diagonal scalar and vector part of the density matrix, respectively. 

The finite-range direct contribution can generally be treated exactly, while the finite-range 
exchange contribution has to be treated using certain approximation. To get similar density 
functional form as that from the zero-range Skyrme interaction, we use in the following 
calculation the density-matrix expansion method [281 129] . By introducing the coordinate 
transformation r = (fx + f*2)/2 and s = n — r*2, the off-diagonal scalar and vector part of 
the density matrix can be expanded in the lower orders as 



s 

p T (r+ -,f- ^) « psL(k T s)p T (r) + g(k T s)s z 



-V 2 p T (f)-T T (r) + -k 2 p T (r) 



(9) 



and 



ssz — * 
Pr(r + - g) ~ 2-?'o(Avs)s x J T (r). (10) 

In the above, k T = (3ir 2 Pt) 1 ^ 3 is the Fermi momentum, and r r and J T are the kinetic density 
and the spin-current density, which can be respectively expressed as 

^(f) = £|V(Mr->)| 2 , (11) 
i , <7 

and 

X(f) = -^5>*,(ry)V0 Ti (r>) x (a'\a\a). (12) 



2 a,a' 



The function psh{k T s) and g(k T s) can be expressed respectively as 

/, \ 3ji(fc T s) 
PSL{k T s) = — — — , (13) 

and 

/, s 35j 3 (k T s) 

9{Ks) = w (14) 

where jx and j'3 are the first- and the third-order spherical Bessel function. 

Applying Eqs. §§§ and (|1U|) to Eq. (0), a Skyrme-like potential energy density functional form 
from the exchange contribution of the Yukawa potential can be obtained with density-dependent 
coefficients as 

H§ (r) = £ {A[p T (r)\ + S^(r)]r T (f) + C[p T {r)] [Vp T (r)] 2 + ^[p r (r)] J 2 (r)}, (15) 

r=n,p 

where the expressions of the coefficients are 

A( Pr ) = -\J P y sL (k T s)V Y (s)d 3 s-^ 2 ) 2 / 3 p 5 T / s B(p T ), (16) 
B(p T ) = - J p T psL(krs)g(k T s)s 2 V Y (s)d 3 s, (17) 

c M - 1^1, , 18 , 

<p{pr) = -~ j ' j 2 (k T s)s 4 V Y (s)ds, (19) 



with V Y (s) = g 2 e~ fls /(4TTs). 

In the standard Skyrme-Hartree-Fock calculation, the finite nuclei properties are obtained by 
solving the following Schrddinger equation 



-V- 



Ti 2 



2m* (r) 



V + U T (r) + W T 



-iV x a) 



(20) 



where e T i is the eigenvalue of the single-particle energy and cj) T i is the nucleon wave function. 
Due to the existence of the Yukawa-type non-Newtonian gravitational potential, the single- 
particle potential U T , the nucleon effective mass m*, and the spin-orbit potential W T are not 
only determined by the Skyrme interaction but modified by the Yukawa potential as well, and 
they can be calculated from the variational principle. The modifications of the single-particle 
potential from the Yukawa potential are from both the direct contribution and the exchange 
contribution 



U? = C/ r YD + C/ r YE , (21) 

vJ° = j *f)$~W=¥[* 1 '> (22) 

uJE = d_A^) + dBM _ dC[pA (VpT? _ 2C(px)v 2 pT + MpAj*. (23) 

We note that the direct contribution is much larger than the exchange contribution. The nucleon 
effective mass and the spin-orbit potential are modified only by the exchange contribution of 
the Yukawa potential, respectively, according to 

h 2 h 2 

2m* 2m* ^ ^ 

and 

H> T Y = 2tp(p T )J T . (25) 



3. Results and discussions 

Using the formulism discussed in Sec. [21 we are now able to study the properties of finite nuclei 
in the presence of the Yukawa-type non-Newtonian gravitational potential. In the following, we 
will discuss the effects from the Yukawa potential by comparing its strength with the Coulomb 
interaction in finite nuclei and then set an upper limit of its strength. 



3.1. Comparing with the Coulomb interaction 

It is instructive to set the strength of the Coulomb interaction as a base line and discuss how 
the finite nuclei properties will be affected if the non-Newtonian potential in finite nuclei is as 
strong as the Coulomb potential. In Fig. [1] we compare the Coulomb potential with the non- 
Newtonian Yukawa potential between two protons with a = —1.24 x 10 36 or —1.24 x 10 34 , and 
A = 1 fm or 10 fm, respectively. We note that a = —1.24 x 10 36 leads to the coupling constant 
<7 2 /47r = 1/137, which is exactly the same strength as the Coulomb potential. Thus, at short 
distances the Yukawa potential is similar to the Coulomb potential while it decreases faster for 
a smaller value of A due to the exponential decay term. The Yukawa potential with only 1% 
strength of that is also shown for reference. 

Next, let's see what would happen in finite nuclei if the Yukawa potential is as strong as the 
Coulomb potential. It is seen in Fig. [2] that the charge density profile of lead nucleus will be 
more diffusive in the presence of the repulsive Yukawa potential, compared to the result from 




Figure 1. Comparing the Coulomb potential with the non-Newtonian Yukawa potential 
between two protons of different strength and length scale parameters at different distances 
of separation. Taken from Ref. |30| . 



MSLO force only. Thus, the repulsive Yukawa potential as strong as the Coulomb potential 
would generally increase the charge radii of finite nuclei, and the effect is larger with a larger A. 
This is quite understandable, and it is seen that the effect becomes negligible if the strength of 
Yukawa potential is 1% of the Coulomb potential. 




Figure 2. Charge density profiles of lead nucleus from Skyrme-Hartree-Fock calculation with 
the MSLO Skymre force only and with an additional non-Newtonian Yukawa potential of different 
strength and length scale parameters. Taken from Ref. [30j . 



3.2. Nuclear Constraint on the Non-Newtonian potential 

From the above discussion, we intuitively know that the Yukawa-type non-Newtonian 
gravitational potential will change the finite nuclei properties unless it is of short range and/or 
much weaker than the Coulomb potential. To avoid modifying the well-determined nucleon- 
nucleon interaction available, we require that the Yukawa potential will not change much the 
finite nuclei properties which have been well described by our framework in nuclear physics, at 
least within the uncertainty from nuclear predictions. Especially, the charge radii and binding 
energies, which can be very accurately measured in experiments, are mostly related to the well- 
constrained isoscalar part of the nucleon-nucleon interaction but are not sensitive to the less 
certain isovector part and/or three-body interaction, and they can now be reproduced by the 
MSLO Skyrme force within about 2%. This 2% uncertainty is thus the largest room for the 
possible existence of the non-Newtonian gravity at femtometer scale. 




Figure 3. The charge radii (r c ) and binding energies per nucleon (B.E.) of 208 Pb, 120 Sn, and 
Ca nuclei as functions of the strength parameter a of the non-Newtonian Yukawa potential 
at length scales A = 1 and 10 fm. The horizontal lines are the mean values of the experimental 
data |3H 132] . a = is the case without non-Newtonian potential, and the green band is the 
largest room for a that the change of the charge radius or binding energy of 208 Pb due to the 
existence of the non-Newtonian potential is within 2%. Taken from Ref. [30J. 

How the charge radii and binding energies from medium to heavy nuclei change with the 
strength parameter a of the Yukawa potential at length scales A = 1 and 10 fm is displayed 
in Fig. In the presence of the repulsive (attractive) Yukawa potential, corresponding to a 
negative (positive) value of a, the charge radii increase (decrease) and the nuclei become less 
(more) bound, which is quite understandable. The change is almost linear to a and more 
sensitive to the strength parameter for A = 10 fm than for A = 1 fm. In addition, the change 
is larger for heavy nuclei 208 Pb due to the finite-range nature of the Yukawa potential. Also 
plotted in the figure is the 2% uncertainty of both charge radius and binding energy of 208 Pb 
from the MSLO force only, which are shown with green bands. The green bands are thus the 
constraints for the strength parameter of the Yukawa potential. 

The constraints on the strength of the non-Newtonian potential at femtometer scale from 
nuclei binding energies and charge radii obtained above are compared with those from other 



log (ji(MeV)) 



40 



35 



CT> 

^ 25 



20 



15 




not allowed 



30 - 




allowed 



(a) n-Pb scattering and X-ray 
from muonic atoms 
(Barbieri era/. 1975) 

(b) n-mirror scattering 
(Nesvizhevsky er a/. 2004) 

(c) spectrum of antiprotonic 

He atoms and n-Pb scattering 
(Pokotilovski 2006) 
I I I I 



(d) n-nucleus forward scattering 
and total cross section 
(Nesvizhevsky er al. 2008) 

(e) np scattering 
(Kamyshkov et al. 2008) 

(f) nuclei charge radii (this work) 

(g) nuclei binding energies (this work) " 



-5 



-10 



-15 



-20 



O 
to 



-15 



-14 



-13 



-12 



-11 



-10 



log (A,(m)) 



Figure 4. Constraints of the non-Newtonian Yukawa potential at femtometer scale from nuclei 
binding energies and charge radii in the present work together with those from neutron scattering 
experiments at larger scales. Taken from Ref. |30] , 



works mostly using neutron scattering experiments [33j [HH [35l [36l E] at larger scales in Fig. [H 
It is found the upper limit is larger at shorter distances, reflecting the increasing difficulty of the 
experiments and the stronger potential at work at smaller scale. In addition, it is found that 
our constraint covers the previously unexplored region and is a smooth expansion of those at 
larger scales. Typically, the constraint from nuclei binding energy can be written as 

log(H) < 1.75/[A(fm)] - 54 + 33.6, (26) 

and that from nuclei charge radius can be written as 

log(|a|) < 1.18/[A(fm)] - 79 + 35.0, (27) 

for A = 1 ~ 10 fm. Since the constraint from the binding energy is stronger, it can be used as 
the nuclear constraint on the non-Newtonian gravity at femtometer scale, and it also leads to 
the upper limit of the boson-nucleon coupling constant in the form of 

log{g 2 ) < 0.10[^(MeV)] - 54 - 3.53, (28) 

if the boson mass [i is between 20 MeV and 200 MeV. 

4. Summary and Outlook 

In this talk we report our recent results of nuclear constraints on non-Newtonian gravity at 
femtometer scale. To do this, we consistently incorporate the Yukawa-type non-Newtonian 
gravitational potential in the Skyrme-Hartree-Fock calculation. The nucleon interaction is 
represented by the well-established MSL0 Skyrme force. By assuming that the Yukawa potential 



will not change the binding energies and charge radii of finite nuclei by 2%, the strength of the 
Yukawa potential is constrained within log(|a|) < 1.75/[A(fm)] ' 54 + 33.6 for A = 1 ~ 10 fm, 
so that the calculated properties of finite nuclei will not be in conflict with the very accurate 
experimental data available. This constraint may serve as a useful reference in constraining 
properties of weakly-coupled gauge bosons and further explorations of possible extra dimensions 
at femtometer scale. 
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